PROCEEDINGS 

OF SCIENCE 



(N 
O 

O 

Q 



Oh. 

CD 



> 

O 

(N 

m 
(N 

(N 



X 



The sphaleron rate at the electroweak crossover 
with 125GeV Higgs mass 



Michela D'Onofrio 

Department of Physics, University of Helsinki and Helsinki Institute of Physics 

P.O. Bo x 64, FI-00014 Helsinki 

E-mail: [niche la . donof rio@helsinki . f i 



Kari Rummukainen 

Department of Physics, University of Helsinki and Helsinki Institute of Physics 

P.O. Box 64, FI-00014 Helsinki 

E-mail: [kari . rummukainen@helsinki . f i 



Anders Tranberg 

Niels Bohr Institute, University of Copenhagen 
Blegdamsvej 17, DK-2100 Copenhagen 
and Helsinki Institute of Physics 
P.O. Box 64, FI-00014 Helsinki 



E-mail: antranbe@nbi.dk 



We measure the sphaleron rate with the physical parameters of the Standard Model. In particu- 
lar, we plug into the calculations the recently found Higgs mass mn = 125 GeV. The sphaleron 
rate tells us about the efficiency of baryon number violation through sphaleron transitions. These 
occur above the electroweak scale ~ 100 GeV and get exponentially suppressed at temperatures 
substantially below the electroweak crossover. The sphaleron rate enters computations of Baryo- 
genesis via Leptogenesis, where it converts non-zero lepton number into non-zero baryon num- 
ber. We simulate the effective electroweak theory on the lattice with multicanonical and real-time 
methods to calculate the sphaleron rate as a function of temperature through the electroweak 
crossover. 



The 30th International Symposium on Lattice Field Theory 
June 24 - 29, 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Lieence. 



http: //pos . sissa . it/ 



The sphaleron rate at the electroweak crossover with 125 GeV Higgs mass 



Michela D'Onofrio 



1. Introduction 

Baryon and lepton numbers are not exactly conserved quantities in the Standard Model, be- 
cause of the axial anomaly, which connects them to the Chern-Simons number of the weak gauge 
field. Vacua in the electroweak theory are labeled by an integer-valued Chem-Simons number 

2 

N CS = J d 3 x f cs = - J d\ e'^Tr (a,F ; * + i^AjA^ . (1.1) 

In the early Universe, for T > 100 GeV, it is possible to move from one vacuum to another sur- 
mounting the potential barrier through sphaleron transitions, thus changing Chern-Simons number, 
B and L according to 

— [B(t)-B(0)]=L i (t)-L i (0)=Ncs(t)-N cs (0), (1.2) 
n G 

with no the number of generations of fermions. At zero temperature, where the potential barrier is 
high, the rate of Chern-Simons number fluctuations, the sphaleron rate, is exponentially suppressed 
and, when the Higgs field expectation value v 3> T, the rate is negligible. In electroweak Baryo- 
genesis scenarios [JXJ] the baryon number of the Universe is generated during the electroweak phase 
transition. However, this scenario does not work in the Standard Model: it requires stronger CP 
violation and first order phase transition, whereas the Standard Model has a smooth crossover 

Nevertheless, the sphaleron rate during the electroweak crossover in the Standard Model is 
relevant for Baryogenesis via Leptogenesis: indeed the sphaleron rate converts lepton asymmetiy 
into baiyon asymmetry. If the lepton asymmetry is generated just before or during the electroweak 
phase transition, how the sphaleron rate shuts off has an effect on the generated baryon number. 

The sphaleron rate has been studied in the broken phase before, but either with unphysical 
Higgs masses f|, or not very deeply in the broken phase Both perturbative calculations 
[§, 0] and lattice simulations |J ||, 10] have been used. 



In this paper we follow the procedure of our previous works [ |1 1| , |12| ] in determining the 
sphaleron rate at the energy range of the electroweak crossover, with the crucial difference that 
now we are able to perform the simulations with the recently found Higgs mass, m# =125 GeV. 
Our results are compared to analytical estimates both in the broken and symmetric phases [§]. 

2. Theory on the lattice 

The thermodynamics of the 4-dimensional electroweak theory is studied in 3 dimensions by 



means of dimensional reduction [13], a perturbative technique that gives the correspondence be- 
tween 4D and 3D parameters. The result is a SU(2) effective theory with the Higgs field and 
gauge field (Fy) 

L = ^FfjFfj + {D4)\D$) +ml^(j> + A 3 (0 t 0) 2 , (2.1) 
and 3D effective parameters g\, A3 and m^. 



Bodeker showed [14] that at leading order in log(l/g) the time evolution of this effective SU(2) 



Higgs model is governed by Langevin dynamics. The latter, though, is very slow on the lattice and 
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can be substituted by any other dissipative procedure, e. g. heat bath. One heat-bath sweep through 
the lattice corresponds to the real-time step [f|] 

At = ^, (2.2) 



with 

-1 3 Ng 2 T 



ln™°+ 3.041 



(2-3) 



An 

where a e \ is the non-abelian color conductivity, which quantifies the current response to infrared 
external fields, N is the dimension of the SU(N) gauge group, and mu is the Debye mass, determin- 
ing the length scale Id ~ 1 /mo ~ 1 /gT. We made use of a 32 3 lattice, with fio = -4- = 9, where g% 
is the 3D gauge coupling and a the lattice spacing. In real-time simulations, for each temperature 
we computed 4 trajectories for every 1000 initial configurations. 

3. Measurement of the sphaleron rate 

In the symmetric phase we make use of canonical Monte Carlo simulations and approach the 
broken phase. At very low temperatures, the rate is highly suppressed and canonical methods 
do not work anymore. We need multicanonical methods, which calculate a weight function that 
compensates the high potential barrier between the vacua, thus allowing transitions. The exact 
value of the sphaleron rate 

rglim ((Afa«)-Afa(0tf) 
t^oo y t 

is obtained, in the broken phase, through a method similar to the one used in [|,|]. 

1. Once done the multicanonical simulations, we obtain the canonical (physical) probability 
distribution of the Chem-Simons number p p h ys .(/Vcs)- 

2. We choose a narrow interval 1/2 — e/2 < Nqs < 1/2 + 6/2 around the point that separates 
vacuum Ncs = from the vacuum jVcs = 1- The relative probability of finding a configura- 
tion here is 

/■ 1/2+6/2 

F(|ATcs-l/2| <e/2)= / dN Pphys (N). (3.2) 

Jl/2-e/2 

This is where we need multicanonical methods, as the probability of being on top of the 
barrier is extremely small, and to get a reliable estimate would take an impractically long 
time with canonical sampling. 

3. Let us now take a random configuration from the canonical distribution but with the con- 
straint 1/2 — e/2 < Ncs < 1/2 + e/2; i.e. near the top of the potential barrier. Starting from 
this configuration, we now generate two real-time trajectories using heat-bath dynamics. 
The trajectories are evolved until the Chem-Simons number falls near a vacuum value. In- 
terpreting one of the trajectories as evolving backwards in time, we can glue the trajectories 
together at the starting point and obtain a vacuum-to-vacuum trajectory. The trajectory can 
either return to the starting vacuum or be a genuine tunneling trajectory. Only the latter-type 
trajectories contribute to the sphaleron rate. 
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4. We can obtain the tunneling rate by measuring |A2Vcs/Af | from the trajectories at the moment 
they cross the value Nqs = 1/2- Here At is the time interval between successive measure- 
ments, and AAfcs the change in Chern-Simons number. This characterizes the probability flux 
throughout the top of the barrier. We obtain the physical time difference from the relation 



between the heat-bath "time" and physical time, equation (2.2) 



5. If the tunneling trajectories would go straight across the top, the ingredients above would 
be sufficient to calculate the total rate. However, typically the trajectories "random walk" 
near the top of the barrier and can cross the value Nqs = 1/2 several times. Because the 
trajectories were chosen starting from a set of configurations near the top of the barrier, this 
leads to overcounting: the evolution could be started at any point the Nqs = 1/2 limit is 
crossed. This can be compensated by calculating a dynamical prefactor 

d = —Y gtunnel , (3.3) 
^traj £y # crossings 

where the sum goes over the ensemble of trajectories, A^j is the number of trajectories, 
^tunnel is if the trajectory does not lead to a change of the vacuum and 1 if it does, and 
(# crossings) is the number of times the trajectory crosses Nqs = 1/2. 

With these ingredients, the sphaleron rate now becomes 

AAfcs 



P(|iVcs-l/2| <e/2) / 



At 



d. (3.4) 



We note that the result is independent of e as long as £ 1. It is also independent of the frequency 
At with which the Chern-Simons number is measured: if we decrease the measurement interval, the 
trajectories become more jagged due to the random-walk nature of the heat-bath updates. This will 
increase the number of the crossings of the value Afcs = 1/2 and hence decrease d. However, the 
latter is completely compensated by a corresponding increase in (|A2Vcs/Af|). If the measurement 
interval At is small enough, random walk arguments imply d oc (A?) 1 / 2 and {\ANcs/At\) ^ (Af) -1 / 2 . 
This is corroborated by the numerical data. Thus, equation (|3.4[) has a well-defined continuum limit. 



4. Results 

We obtain the sphaleron rate T/T A and the expectation value for the Higgs field (0 2 ) for Higgs 
mass niH =125 GeV. Figure [I] shows the Higgs field expectation value behavior as a function of 
temperature. We start from the "symmetric phase" with canonical Monte Carlo simulations and 
lower the temperature to reach the "broken phase", where we switch to multicanonical simulations. 
We can see the Higgs field assuming a non-zero value when approaching the broken phase, and the 
transition from the two methods occurring smoothly. 

The sphaleron rate as a function of temperature is shown in Figure ||. Here again we perform the 
simulations with canonical Monte Carlo at high temperatures and continue with multicanonical 
methods when reaching the cold broken phase. The sphaleron rate changes from its asymptotic 
value to become exponentially suppressed at very low temperatures. The canonical and multi- 
canonical methods are in good agreement. The theoretical curves, used to compared our results, 
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Figure 1: The Higgs expectation value (0 2 ) for Higgs mass of 125 GeV as a function of temperature. The 
high-temperature canonical and low-temperature multicanonical results match beautifully in the transition 
region. 



were obtained separately for the broken and symmetric phases, through perturbative calculations 
in®. 



5. Conclusion 



We improved the previous estimates for the sphaleron rate and determined its behavior from 
the symmetric to the broken phase, through the electroweak crossover. Our results are in agreement 
with previous estimates in the symmetric phase, and in the broken phase the slope of our curve is 
the same as in the analytic calculation ®. 

Even though the Standard Model has a too weak source of CP-violation in the quark sector, 
Baryogenesis might still be viable through lepton number violating processes. The sphaleron rate 
plays an important role in Leptogenesis, as the conversion of lepton to baryon number depends on 
it, and it is therefore important to know its size rather accurately. 
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Figure 2: The sphaleron rate for a Higgs mass of 125 GeV. The high-temperature canonical and low- 
temperature multicanonical results again match very well in the transition region. Also shown are previous 
high-temperature estimates ~ 5 x 10~ 7 T 4 , (top, horizontal line) and perturbative calculations in the low- 
temperature phase (bottom, wide band), both from [||]. 
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